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Abstract 



A new class of graphical models capturing the dependence structure of events 
that occur in time is proposed. The graphs represent so-called local independencies, 
meaning that the intensities of certain types of events are independent of some (but 
not necessarilly all) events in the past. This dynamic concept of independence is 
^ ^ , asymmetric, similar to Granger non-causality, so that the corresponding local inde- 

C/^ ■ pendence graphs differ considerably from classical graphical models. Hence a new 

• notion of graph separation, called 5-separation, is introduced and implications for 

^ ' the underlying model as well as for likelihood inference are explored. Benefits re- 

garding facilitation of reasoning about and understanding of dynamic dependencies 
as well as computational simplifications are discussed. 
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O '. 1 Introduction 



I Marked point processes are commonly used to model event history data, a term originating 
from sociology where it is often of interest to investigate the dynamics behind events 
such as finishing college, finding a job, getting married, starting a family, durations of 
unemployment, illness etc. But comparable data situations also occur in other contexts, 
e.g. in survival analyses with intermediate events such as the onset of a side effect, change 
of medication etc. (cf. Keiding, 1999). Longitudinal studies and the careful analysis of the 
underlying processes are crucial for gaining insight into the driving forces of inherently 
dynamic systems, but having to deal with the multidimensionality as well as with the 
dynamic nature of these systems makes this a very complex undertaking. 

Graphical models deal with complex data structures that arise whenever the interre- 
lationship of variables in a multivariate setting is investigated. Over the last two decades, 
they have proven to be a valuable tool for probabilistic modelling and multivariate data 
analysis in such different fields as expert systems and artificial intelligence (Pearl, 1988; 
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Cowell et al., 1999; Jordan, 1999), hierarchical Bayesian modelhng (cf. BUGS project 



http://www.mrc-bsu.cain.ac.uk/bugs/welcome.shtml), causal reasoning (Pearl, 2000; 



Spirtes et al., 2000), as well as sociological, bio-medical, and econometric applications. For 
overviews and many different applications see for instance the monographs by Whittaker 
(1990), Cox and Wermuth (1996), Edwards (2000). 

While the 'classical' graphical models are concerned with representing conditional 
independence structures among random variables, variations have been proposed to deal 
with feedback systems (Sprites, 1995; Koster, 1996) but are still based on cross-sectional 
data. The application of graphical models to truly time-dependent data, such as event 
histories or time series is only slowly getting on its way. Dahlhaus (2000) has proposed 
graphical models for multivariate time series, but his approach does not capture how the 
present or future of the system depends on or is affected by the past. Instead, Eichler 
(1999, 2000) uses graphs to represent Granger-causality, which is a dynamic concept of 
dependence. For continuous time one approach, called Dynamic Bayesian Networks, is to 
discretise time and provide directed acyclic graphs that encode the independence structure 
for the transitions from ttot + 1 (Dean and Kanazawa, 1989). The approach proposed by 
Nodelman et al. (2002, 2003) comes closest to the kind of graphs that we will consider. 
In their Continuous Time Bayesian Networks they represent a multistate Markov process 
with nodes corresponding to subprocesses and edges corresponding to dependencies of 
transition rates on states of other subprocesses. 

In this paper, we propose and investigate the properties of graphs that represent 
so-called 'local independence' structures in event history data. The basic idea of local 
independence is that, once we know about specific past events, the intensity of a consid- 
ered future event is independent of other past events. It has been developed by Schweder 
(1970) for the case of Markov processes and applied e.g. in Aalen et al. (1980). A gen- 
eralisation to processes with a Doob-Meyer decomposition can be found in Aalen (1987) 
who focusses on the bivariate case, i.e. local dependence between two processes. Here 
we extend this approach to more than two processes. The analogy of the bivariate case 
to Granger-non-causality has been pointed out by Florens and Fougere (1996), see also 
Comte and Renault (1996). Note that the notion of 'local independence' used by AUard 
et al. (2001) is a different one. 

We first set out the necessary notations and assumptions for marked point processes 
in section 12.11 followed by the formal definition of local independence in section 12.21 the 
emphasis being on the generalisation to a version that allows to condition on the past 
of other processes and hence describes dynamic dependencies for multivariate processes. 
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Section 13.11 defines graphs that are appropriate to represent local (in)dependence. The 
main results are given in section 13. 2[ The properties of local independence graphs are 
investigated. In particular we prove that a new notion of graph separation, called 6- 
separation, can inform us about independencies that are preserved after marginalising 
over some of the processes. In section 13.31 it is shown how the likelihood of a process with 
given local independence graph factorises and implications are discussed. The potential of 
local independence graphs is discussed in section H] and proofs are given in the appendix. 

2 Local independence for marked point processes 

Marked point processes are briefly reviewed in section 12.11 using the notation of Andersen 
et al. (1993). In section the concept of local independence is explained in some detail. 

2.1 Marked point processes and counting processes 

Let S = {ei, . . . jCk}, K < oo, denote the (finite) mark space, i.e. the set containing all 
types of events of interest for one observational unit, and T the time space in which the 
observations take place. We assume that time is measured continuously so that we have 
T = [0, r] or T = [0, r) where r < oo. The marked point process (MPP) Y consists of 
events given by pairs of variables (T^, Es), s = 1, 2, . . . , on a probability space (fi, JF, P) 
where G T, < Ti < T2 . . . , are the times of occurrences of the respective types of 
events Eg G £. Assume that the MPP is non-explosive, i.e. only a finite number of events 
occurs in the time span T. The mark specific counting processes Nk{t) associated with an 
MPP are then given by 

Nkit) = J2HEs = ek}, k = l,...,K. 

Ts<t 

We write N = (A''i, . . . , Nk) for the multivariate counting process, and N^, Ac {1, . . . , K}, 
for the vector {Nk)keA, calling a subprocess, with Ny = N. 

To investigate dependencies of the present on the past it will be important to have 
some notation for the history of some subset or all of the processes involved. Hence we 
denote the internal filtration of a marked point process by J-'t = (^{{Tg, Eg) \ Tg < t, Eg C 
S} which is equal to a{{Ni{s), . . . , Nk{s)) \ s < t}, while for A C {1, . . . , K} we define 
the filtrations of a subprocess as J-'f^ = (t{Na(s) | s < t}; in particular JF^^ is the internal 
filtration of an individual counting process N^. 

Under rather general assumptions (cf. Fleming and Harrington, 1991, p. 61), a Doob- 
Meyer decomposition of Nk{t) into a compensator and a martingale exists. Both these pro- 
cesses depend on the considered filtration which is here taken to be the internal filtration 
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of the whole MPP Y. We will assume throughout that all the jF^-compensators are ab- 
solutely continuous and predictable so that intensity processes Xk(t) exist, which are taken 
to be predictable versions of the derivatives of the compensators, i.e. Afc(t) = \k{s)ds. 
Heuristically we have (Andersen et al., 1993, p. 52) 

Xk{t)dt = E{Nk{dt)\Tt-)- (1) 

More formally this means that the differences A^^^ — A^ are jF^-martingales. 

An interpretation of property ([1]) is that given the information on the history of the 
whole MPP up to just before time t, Xk{t)dt is our best prediction of the immediately 
following behaviour of N^. Note that for the theory developed in this paper the setting 
can slightly be generalised, not requiring absolute continuity of compensators (cf. Didelez, 
2000), as might be relevant when certain types of events can only occur at fixed times. 

It will be important to distinguish between the ^(-intensity based on the past of the 
whole MPP, and the jF^^^-intensities based on the past of the subprocess on marks in A. 
The latter can be computed using the innovation theorem (Bremaud, 1981, pp. 83), and 
a way of doing so, especially relevant to our setting, is given in Arjas et al. (1992). 

The following is a standard assumption in counting process theory but we want to 
highlight it as it plays a particularly important role for local independence graphs. 

Assumption 2.1 No jumps at the same time 

The jFf-martingales — J Xk{s)ds are assumed to be orthogonal for k G {1, . . . ,K}, 
meaning that none of A'^i, . . . , Nk jump at the same time. This is implied by the above 
assumption that all compensators are absolutely continuous if in addition no two counting 
processes Nj and A^^ are counting the same type of event. 

Assumption 12. II might be violated, e.g. when investigating the survival times of cou- 
ples and there is a small but non-zero chance that they die at the same time, in a car 
accident for instance. The reason for imposing this assumption is that we want to explain 
dependencies between events by the past not by common innovations. If one wants to 
allow events to occur at the same time, then such a simultaneous occurrence defines a 
new mark in the mark space. 

Note that general multi-state processes can be represented as marked point process 
with every transition between two states being a mark. This is explored in more detail in 
Didelez (2007) for Markov processes. 

2.2 Local independence 

The bivariate case is defined as follows (cf. Aalen, 1987). 
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Definition 2.2 Local independence (bivariate). 

Let Y be an MPP with S = {ei, 62} and A'^i and the associated counting processes on 
{Q,J-',P). Then, A'^i is said to be locally independent of N2 over T if Xi{t) is measurable 
w.r.t J-'l for all t E T. Otherwise we speak of local dependence. 

The process A^^i being locally independent of N2 is symbolised by 7^ ^i- Inter- 
changeably we will sometimes say that ei is locally independent of 62, or 62 7^61. 

The essence of the above definition is that the intensity Ai(t), i.e. our 'short-term' 
prediction of A'^i, remains the same under the reduced filtration J-'^ as compared to the 
full one J^f This implies that we do not lose any essential information by ignoring how 
often and when event 62 has occurred before t. One could say that if 7^ ^1 then the 
presence of A''i is conditionally independent of the past of N2 given the past of A''i, or 
heuristicallj|3 

iVi(t)^^M^,l, (2) 

where AAL B\C means 'A is conditionally independent of B given C" (cf. Dawid, 1979). 
For general processes, this is a stronger property than local independence but it holds for 
marked point processes with Assumption 12.11 as their distributions are determined by the 
intensities. Note that ([2]) does not imply that for m > 0: A^i(t + M)_LLJF?_|JF^L — hence the 
name local independence. Also, Ni(t)ALN2{t) will only hold if the two processes are mu- 
tually locally independent of each other. Without Assumption 12.11 the jF^^-measurability 
of Xi(t) would, for instance, trivially be true if ei = 62 but, in such a case, we would not 
want to speak of independence of ei and 62- 

Example 1. Skin disease. 

In a study with women of a certain age, Aalen et al. (1980) model two events in the life 
of an individual women: occurrence of a particular skin disease and onset of menopause. 
Their analysis reveals that the intensity for developing this skin disease is greater once the 
menopause has started than before. In contrast, and as one would expect, the intensity for 
onset of menopause does not depend on whether the person has earlier developed this skin 
disease. We can therefore say that menopause is locally independent of this skin disease 
but not vice versa. Note that, in whatever way the onset of skin disease and menopause 
are measured, it is assumed that they do not start systematically at exactly the same 

^Notc that (I2I) is an informal way of saying that Ni{dt) is conditionally independent of {Ts\Ts < 
t, Eg = 2,5 = 1,2,...} given {Ts\Ts < t, i?s = 1, s = 1, 2, . . .}. This and similar statements later, like ([3|), 
(fT2l) or (fT5|) . should be interpreted correspondingly. 
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time, corresponding to the above 'no jumps at the same time' assumption. 

Let us now turn to the case of more than two types of events. This requires condi- 
tioning on the past of other processes as follows. 

Definition 2.3 Local independence (multivariate). 

Let N = {Ni, . . . , Nk) be a multivariate counting processes associated with an MPP. 
Let further A, B, C be disjoint subsets of {1, ... , K}. We then say that a subprocess 
is locally independent of given 'Nc over T if all jFj'^^'^'-"^-intensities A^, k E B, are 
measurable with respect to J^f^^' for all t E T. This is denoted by 7^ | Nc or 
briefly A -/^ B \ C . Otherwise, is locally dependent on given Nc, i.e. A ^ B \ C . 
If C = then B is marginally locally (in)dependent of A. 

Conditioning on a subset C thus means that we retain the information about whether 
and when events with marks in C have occurred in the past when considering the intensi- 
ties for marks in B. If these intensities are independent of the information on whether and 
when events with marks in A have occurred we have conditional local independence. In 
analogy to ([2]) the above definition of multivariate local independence is with Assumption 
12. II equivalent to 

NBit)ALJ^f^\J^t^y^ yteT. (3) 

Example 2. Home visits 

This example is not taken from the literature but is inspired by real studies (e.g. Vass 
et al., 2002, 2004). In some countries programmes exist to assist the elderly through reg- 
ular home visits by a nurse. This is meant to reduce unnecessary hospitalisations while 
increasing the quality of life for the person. It is hoped that such a programme increases 
the survival time. The times of the visits as well as the times and durations of hospitali- 
sation are monitored. In addition, it is plausible that the underlying health status of the 
elderly person may also affect the rate of hospitalisation and predict survival. This inter- 
play of events for an individual elderly person can be represented as an MPP if 'health 
status' is regarded as a multistate process. Assume that the timing of the home visits is 
determined externally, e.g. by the availability of nurses which has nothing to do with the 
patient's development, i.e. the visits are assumed locally independent of all the remain- 
ing processes. It might then be of interest to investigate whether the visits affect only 
the rate of hospitalisation directly, i.e. whether survival is locally independent of the vis- 
its process given the hospitalisation and health history or even given only a subset thereof. 
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As can easily be checked, local (in) dependence needs to be neither symmetric, re- 
flexive nor transitive. However, since in most practical situations a subprocess depends 
at least on its own past we will assume throughout that local dependence is reflexive. An 
example for a subprocess that depends only on the history of a different subprocess and 
not on its own history is given in Cox and Isham (1980, p. 122). 

In order to see the relation with local independence, we briefly review Granger non- 
causality (Granger, 1969). Let Xy = {Xyit) \ t e Z} with Xy(t) = {Xi{t), . . . , Xxit)) 
be a multivariate time series, where V = {1, . . . , K} is the index set. For any A G V we 
deflne Xa = {XA{t)} as the multivariate subprocess with components Xa, a E A. Further 
let XA{t) = {Xa{s) \ s < t}. Then, for disjoint subsets A, B C V we say that Xa is 
strongly Granger-non causal for Xb if 

XBit)ALXAit-l) \Xv\Ait-l), 

for all t E Z. The interpretation is similar as for local independence, i.e. the present value 
of Xb is independent of the past of Xa given its own past and the one of all other com- 
ponents C = V\{A U B), in analogy to Also note that the above does not imply that 
Xsit + u)ALXA{t — 1) I Xv\A{t — 1) for u > 0, again analogous to local independence. 
Eichler (1999, 2000) investigates a graphical representation and rules to determine when 
the condition XY\A(t — 1) can be reduced to proper subsets Xc(t — 1), C C V\A. 

Finally, let us indicate how the deflnition of local independence can be generalised 
to stopped processes. This is relevant when there are absorbing states such as death. In 
that case all other events will be locally dependent on this one because all intensities 
are zero once death has occurred. However, the dependence is 'trivial' and not of much 
interest. Let T be a J?-t-stopping time and let N"^ = {N^, . . . ,N]^) he the multivariate 
counting process stopped at T. Then the intensities of N'^ are given by X]^, k E V, and 
local independence can be formalised as follows. 

Definition 2.4 Local independence for stopped processes. 

Let N"^ = {NJ' , . . . , N]^) he a multivariate counting processes associated to an MPP and 
stopped at time T. Then we say that A -f^B \ C if there exist jF^^'-"^-measurable processes 
Afc, ke B, such that the J^^^^-^^^— intensities of are given by Xl{t) = Xk{t)l{t < T}, 
keB. 

The local independencies in a stopped process have to be interpreted as being valid 
as long as t < T. 
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3 Local independence graphs 



We first give the definition of local independence graphs and then investigate what can 
be read off these graphs. 

3.1 Definition of local independence graphs 

An obvious way of representing the local independence structure of an MPP by a graph 
is to depict the marks as vertices and to use an arrow as symbol for local dependence as 
in the following small example. 

Example 1 ctd. Skin disease. 

The local independence graph for the relation between menopause and skin disease is very 
simple, cf . Figure [1] Notice that even with this simple example there is no way of express- 
ing the local independence using a classical graph based on conditional independence for 
the two times Ti = 'time of occurrence of skin disease' and T2 = 'time of occurrence of 
menopause' as these are simply dependent. 

menopause skin disease 

• >m 

Figure 1: Local independence graph for skin disease example. 

For general local independence structures, we will have directed graphs that may 
have more than one directed edge between a pair of vertices, in case of mutual local de- 
pendence, and that may have cycles. More formally, a graph is a pair G = (y,E), where 
V = {1, . . . , K} is a finite set of vertices and is a set edges. The graph is said to be 
directed if E G {(j, k)\j, k G V,j ^ k}. Later we will also need the notion of an undirected 
graph where E C {{j,k}\j,k EV,jj^ k}. Undirected edges {j,k} are depicted by lines, 
j — k, and directed edges (j, k) by arrows, j — > k. If (j, k) E E and {k,j) G E this is 
shown by stacked arrows, j < — ^ k. 

The following property (jlj) is called the pairwise dynamic Markov property, where we 
say dynamic to emphasise the difference to graphs based on conditional independence. 

Definition 3.1 Local independence graph 

Let Ny = (A'^i, . . . , Nx) be a multivariate counting process associated to an MPP Y with 
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mark space E = {ei, . . . , ex}- Let further G = (V, E) be a directed graph, V = {1, . . . , K}. 
Then, G is called a local independence graph of Y if 



for all j,keV: (j, fc) ^ E {j} ^{k}\V\{j, k}. 



(4) 



Home 
visits 



Hospitalisation 



Home 

visits Hospitalisation 



»-1-->^ 

leath Health status 




*Z1 

•< 4 



Death Health status 



(a) 



(b) 



Figure 2: Home visits example; (a) for the whole process; (b) for the stopped process. 
Example 2 ctd. Home visits 

The graph in Figure [2](a) is for the whole process, while (b) shows the local independencies 
for the stopped process (stopping when death occurs). There are no arrows into 'Home 
visits' representing that the rate of visits is locally independent of 'Hospitalisation' given 
'Health status' as well as of 'Health status' given 'Hospitalisation' (while the person is still 
alive), reflecting that the visits are determined externally. The latter local independence 
might be violated if the nurses, on their own account, increase the frequency of their visits 
when they notice that the person's health is deteriorating. The graph further represents 
that survival is locally independent of the visits given hospitalisation and health history, 
and that the health process is also locally independent of the visits given hospitalisation 
history (while the person is still alive obviously). These absent edges could reflect the 
null-hypothesis when investigating whether the visits affect survival in other ways than 
through changing the rate of hospitalisation. 

3.2 Dynamic Markov properties 

The local independence graphs as defined above allow (under mild assumptions) more 
properties to be read off, concerning the dependence structure, than just those given by 
dl]). We may in particular query the graph with the aim of dimension reduction, i.e. with 
questions about which other processes can be ignored while investigating certain local 
independencies. The local dynamic Markov property addressed in section 13.2.11 tells us 
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about the immediately relevant information when considering a single mark and corre- 
sponding Nk. Further, the global dynamic Markov property in section [3.2.21 gives graphical 
rules to identify when the separating set itself can be reduced, i.e. when in (jl]) we do not 
need to condition on all V\{j, k} but just on a true subset. For this we need the notion 
of (5-separation also introduced in section I3.2.2[ 

Some more graph notation will be required. A path between two nodes is defined in 
the obvious way (the formal definition is given in Appendix IA.2I) : we distinguish between 
undirected paths for undirected graphs, directed paths, preserving the direction of edges, 
for directed graphs and trails for connections in directed graphs that do not preserve the 
direction. For directed graphs we further require the following almost self-explanatory 
notations. If a — > b then a is called a parent of b and 6 is a child of a (if a < — ^ b then a 
is both, a child and a parent of b); pa(y4) denotes the set of all parents of nodes in A G V 
without A itself, and ch(74) analogously the set of children of A. The set cl(y4) = pa(y4)U74 
is called the closure of A. If there is a directed path from a to 6 then a is an ancestor of 
b and 6 is a descendant of a; the corresponding set notation is a.n{A) and de{A) (always 
excluding A itself). Consequently, nd(A) = V\{de{A) U A) are the non-descendants of A. 
If pa(yl) = 0, then A is called ancestral. In general, An(yl) is the smallest ancestral set 
containing A, given by A U an{A). 

3.2.1 Local dynamic Markov property 

Definition 3.2 Local dynamic Markov property 

Let G = (V, E) be a directed graph. For an MPP Y the property 

for aWkeV : V\c\{k) rh{k} \ pa(A;), (5) 

is called the local dynamic Markov property w.r.t. G. 

In other words, property 1^ says that every jF^-intensity is J^-f^^'^^-measurable, 
which clearly implies that for any ancestral set A the intensity Xa is jF/^-measurable. 
This property could for instance be violated if two components in pa,{k) were a.s. identical 
which is however prevented by the orthogonality assumption 12.11 As shown in Appendix 
\A.3\ the exact condition for the property (jS]) to follow from (jl]) is that 

jr^ n J^f = J^f^^ V A, 5 C 1/, Vt G T, (6) 

where we define JF® = {0,fi}. Property is called 'conditional measurable separability' 
(Florens et al., 1990) and formalises the intuitive notion that the components of N are 
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'different' enough to ensure that common events are necessarily due to common compo- 
nents. 

Example 2 ctd. Home visits 

Let us consider the question whether the four processes are sufficiently different to ensure 
dH]). If the health process is measured in a way such that it is determined by the number 
and duration of past hospitalisation, not taking any other information into account, this 
assumption might be violated. However, it makes sense and we will assume for this ex- 
ample that the 'Health status' reflects more aspects of a person's health than just past 
hospitalisations. Then it seems plausible that ([6]) is satisfied as the other processes are 
clearly capturing different information anyway. Consequently we can use the local dy- 
namic Markov property to read off that the visits process is locally independent of both, 
hospitalisation and health status (while the person is still alive). 

3.2.2 (5— separation and the global dynamic Markov property 

In undirected graphs we say that subsets A, B C V are separated by C C \^ if any path 
between elements in A and elements in B is intersected by C. This is symbolised by 
AALgB\C. In classical graphical models every such separation induces conditional inde- 
pendence between A and B given C regardless of whether {A, B, C) is a partition of V 
or not. This can obviously lead to considerable dimension reduction if C is chosen mini- 
mally and the graph is sparse. To obtain a similar result for local independence graphs we 
require a suitable notion of separation called (5-separation, introduced below after some 
more graph notation. 

The moral graph is given by inserting undirected edges between any two vertices 
that have a common child (if they are not already joined) and then making all edges 
undirected (two directed edges between a pair of nodes are replaced by one undirected 
edge). This procedure of moralisation will also be applied to an induced subgraph Ga, 
A (Z V, defined as {A, Ea) with Ea the subset of E containing only edges between pairs 
of nodes in A. Finally, for i? C V, let denote the graph obtained by deleting all 
directed edges of G starting in B. 

Definition 3.3 6 -Separation. 

Let G = (y, E) be a directed graph. Then, we say for pairwise disjoint subsets A, B,C G V 
that G 6-separates A from B in G if AALg B\G in the undirected graph (G^Aii{Aubuc))'" 
(the case of non-disjoint A, B, G is given in Appendix IA.2p . 
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Note that except for the fact that we delete edges starting in B, which makes S- 
separation asymmetric, the definition parallels the one for DAGs. This initial edge deletion 
can heuristically be explained by the fact that we want to separate the present of B from 
the past of A and hence we disregard the 'future' of B which is where the edges out of B 
point to; for the same reason only the ancestral set An{A U -B U C) is considered. As for 
DAGs the insertion of moral edges is necessary whenever we condition on a common 'child' 
due to a 'selection effect' by which two marginally independent variables (or processes) 
that affect a third variable (or process) become dependent when conditioning on this third 
variable. Further properties of 5-separation are discussed in Didelez (2006). 

Definition 3.4 Global dynamic Markov property 

Let Ny = (A''i, . . . , Nk) be a multivariate counting process associated to an MPP Y and 
G = {V, E) a directed graph. The property that 

for all disjoint A,B,G : G 5-separates A from B m. G =^ A -{^B \ G. (7) 

is called the global dynamic Markov property w.r.t. G. 

The significance of the global Markov property is that it provides a way to verify 
whether a subset G C V\{A U B) is given such that A -/^B\G, i.e. local independence is 
preserved even when ignoring information on the past of processes in V\{A U B U G). Of 
course, ([7]) is only meaningful if it can be linked to the definition of local independence 
graphs addressed next. 

Theorem 3.5 Equivalence of dynamic Markov properties 

Let y be a marked point process and G = (V, E) a directed graph. Under the assumption 
of and further regularity conditions (cf. Appendix IA.3p . the pairwise, local and global 
dynamic Markov properties, i.e. (jl]), and ([7j) are equivalent. 

The proof is given in Appendix IA.3[ 
Example 2 ctd. Home visits. 

The underlying health status of an elderly person may be difficult to measure accurately in 
practice. Let us therefore investigate the local independence structure when ignoring this 
underlying process altogether, in particular consider the question of whether from Figure 
[2] we can infer that survival is still locally independent of 'Home visits' given only the 
hospitalisation but ignoring the health process. Graphically this means we have to check 
whether the node 'Hospitalisation' alone separates 'Home visits' from 'Death'. As can be 
seen from the corresponding moral graph (for the stopped process) in Figure [3] this is not 
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the case. Hence, even though the home visits are assumed to be determined externally in 
Figure [2] and do not affect survival directly, ignoring the underlying health process may 
lead to a 'spurious' local dependence of survival on the home visits. The reason is that, 
for instance, a history of hospitalisation with preceding home visit predicts survival differ- 
ently from a hospitalisation without preceding home visit — the former might mean that 
the health was especially bad and hence hospitalisation was necessary, while the latter al- 
lows minor health problems that could have been treated by a nurse who was not available. 

Home 

visits Hospitalisation 




Death Health status 
Figure 3: Moral graph for home visits example. 

Notice that intuitively it is clear that if the intensity of the visits depended on the 
underlying health status, i.e. if there was a directed edge from 'Health' to 'Visits', we 
could talk of confounding. Hence it is rather surprising that even when the frequency of 
the visits is controlled externally we may find a spurious dependence. For the discrete 
time case Robins (1986, 1997) has demonstrated that nevertheless in a situation like Fig- 
ure [2] we can draw causal conclusions even when no information on the underlying health 
process is available. However, standard methods that just model the intensity for survival 
with time varying covariates for the times of previous home visits and hospitalisations 
will typically give misleading results due to the conditional association between 'Home 
visits' and 'Death' given 'Hospitalisation'. 

Example 3. Chemo therapy cycles 

To see a more complex example of a local independence graph consider a hypothetical 
study (which is inspired by real studies) where early stage breast cancer patients are 
observed over a period of several months during which they receive at least one but 
usually more cycles of chemo therapy. The tumour size is monitored through palpation. 
The doctors will consider removing the tumour by surgery if the size does not decrease and 
surgery is almost certain if it increases. Furthermore, the chemo therapy may be delayed 
or discontinued if the patient shows a toxic reaction following the treatment but also if the 
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Toxic Chemo 
reaction therapy Surgery 




Anxiety Tumour Death 
size 

Figure 4: Local independence graph for chemo therapy example (stopped process). 



patient requests a delay, possibly due to an increased state of anxiety. Except for tumour 
size all processes count one type of event that can occur once or more often. Tumour size 
is measured categorically depending on the number and palpable size of lesions and can 
be regarded as a multistate process. Figure S] shows a hypothetical local independence 
structure. For instance it assumes that survival locally depends on the tumour size and 
whether surgery has taken place, but once this information is given none of the other 
processes are relevant for the intensity of death. Note that this particular assumption 
could plausibly be violated because toxic reactions and anxiety may reflect other health 
problems, but for simplicity we will assume that all patients are 'healthy' except for the 
breast cancer so that this violation is excluded. 

Figure O shows the different moral graphs constructed from Figure H] to investigate 
5-separations. Graph (a) shows (5-separation from the node 'Death' allowing us to read 
off, for instance, that 'Chemo therapy' is not (5-separated from 'Death' by 'Tumour size' 
alone reflecting that chemo therapy predicts survival if surgery history is ignored. This 
is plausible because knowing that for example a decrease in tumor size was preceded by 
a treatment cycle is informative for surgery, making it less likely than without preced- 
ing chemo therapy; and whether surgery has taken place, in turn, predicts the survival 
chances. As 'Anxiety' is problematic to observe and measure we may further be interested 
in the question of when it can be ignored. We see from the graph (a) that 'Death' is lo- 
cally independent of 'Anxiety' given either the set {'Surgery', 'Tumour size'} or {'Chemo 
therapy', 'Tumour size'}, the latter implying that once we know the chemo therapy his- 
tory in addition to the development of the tumour size then anxiety will not inform us 
any further about the intensity for death regardless of whether surgery and toxic reaction 
history is known or not. But note that even though 'Anxiety' does not affect 'Tumour size' 
directly the latter has to be part of the separating set. Figure |5l^b) shows that 'Anxiety' 
is (5-separated from 'Surgery' by any set that includes 'Tumour size' and similarly for 
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Figure 5: Different moral graphs for chemo therapy example. Dotted edges have been 
added due to a 'common child' in Figure H] and circles indicate that arrows out of these 
nodes have been deleted before moralising. 

(c) that it is 5-separated from 'Toxic reaction' by any set that includes 'Chemo therapy' 
— in these two cases (5-separation does not tell us more than the local dynamic Markov 
property (jS]). From graph (d) we see that 'Anxiety' itself is locally independent of 'Chemo 
therapy' and 'Tumour size' given 'Toxic reaction' and of 'Surgery' given either 'Toxic 
reaction' or the set {'Chemo therapy', 'Tumour size'}. 

3.3 Likelihood factorisation and implications 

In order to discuss properties and implications for the likelihood for graphical MPPs we 
will regard the data consisting of times and types of events (ti, ei), (^2, 62), ... , (tn, e„) as 
a realisation of the history process Ht = {(Tg, Es)\Ts < t}. As for filtrations, Hf- denotes 
the strict pre-t history process. Additionally, Hf, Ac {1, . . . , K}, defined as 

Hf = {{Ts, Es) \ Ts < t and 3 k e A : Es = Ck, s = 1,2 .. .} 

denotes the history process restricted to the marks in A. Any set of marked points for 
which it holds that tg = tu,s ^ u, implies Cg = can be a history, i.e. a realization of 
Ht. Note that (up to completion by null sets) the different filtrations can be regarded as 
being generated by the history processes, i.e. J^^ = a{Hf^}, A G {1, . . . , K}. 
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Before deriving the likelihood for a given local independence graph, we recall it for 
the general case. Based on the mark specific intensity processes Xk{t) the corresponding 
crude intensity process is given by X(t) = Ylk=i ^k(t)- This is the intensity process of the 
cumulative counting process J2k-^i^- '^^^ likelihood process L{t\Ht) is then given as 

L{t\Ht) = W XeXTs) ■ exp (^-J\{s)ds^ . (8) 

To see how the likelihood is affected by G being a local independence graph of Y, we first 
rewrite ([8]) as follows: 

k=lTs<t \ k=l J 



n n ^k{Ts(k)) ■ exp f - / 



\k(s) ds 



where Ts(^k) with -E's(fc) = ^k are the occurrence times of mark Cfc. The inner product of the 
above can be regarded as the mark specific likelihood and is denoted by Lk{t\Ht). Now, by 
the definition of a local independence graph and the equivalence of the pairwise and local 
dynamic Markov properties under condition ([6]) we have that Afc(s) is J^-f^'-'^^-measurable, 
where c\{k) is the closure of node k. Hence, it follows that 

Lk{t\Ht) = Lkmf-'^), (9) 

i.e. the mark specific likelihood Lk based on the whole past remains the same if the 
available information is restricted to how often and when those marks that are parents of 
Cfc in the graph and Ck itself have occurred in the past, symbolised by Hf^^^\ 
It follows that under (EI) the likelihood factorises as 

rCl(fc)^ 

LJk\t^\-Lj 

kev 

which parallels the factorisation for DAGs where the joint density is decomposed into the 
univariate conditional distributions given the parents. Here, we replace the parents by the 
closure because we also need to condition on the past of a component itself which is not 
required in the static case. 

Example 2 ctd. Home visits 

From Figure [2](b) we obtain the following factorisation 



L{t\H,) = l[Lk{t\Hf^'^), (10) 
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for t < time of death, where vi, ho, hs, d stand for 'visits', 'hospitahsation', 'health status' 
and 'death' respectively. 

Two consequences of the above factorisation regarding the relation of local and con- 
ditional independence are given next. 

Theorem 3.6 Conditional independencies 

For an MPP with local independence graph G and disjoint A,B,C C V, such that C 
separates A and B, i.e. AALgB \ C, in (G'An(AuBuc))™! have 

J^fALJ^flJ^^ \fteT. (11) 

The proof is given in appendix IA.4[ The graph separation that A,B,C have to satisfy 
for ([TT]) implies that for each k e C the J'/^^^^'^-intensity is either Tf^'^- or Tf^'^- 
measurable, otherwise C could not separate A and B in the moral graph. Also, of course, 
we have that for each k ^ A (resp. 5) the jF^"^^^'-"-^-intensity is T^^'" (resp. Tf^'"^- 
measurable; otherwise there would be edges linking A and B in the graph and they could 
not be separated. A property similar to (llip has been noted by Schweder (1970, Theorems 
3 and 4) for Markov processes. With a similar argument we can reformulate ([2]) and ([3]): 
for any 5 C we have 

NB(t)X-7v^'^^^^ I Tf}^\ (12) 
i.e. the present of is independent of the past of Ny^^^^^) given the past of N,,^^). 

Example 3 ctd. Chemo therapy cycles 

Let A =' Surgery', B =' Toxic reaction' and C = {'Chemo therapy', 'Tumour size'}. Then 
(^An(AuBuc))'" is the same as Figure [5l^a), where the node 'Death' could be omitted as 
we are conditioning on the patient being alive anyway, and indeed A and B are separated 
by C. With fITT]) we can infer that at any time t (before death) the whole surgery history, 
i.e. whether and when surgery has taken place before t, is independent of whether and 
when toxic reactions have occured given we know the tumour size development up to t 
and when chemo therapy has been administered. Note that as mentioned above we have 
for the nodes in C that the jF^'^'-^-^^'^-intensity for 'Tumour size' is T^^'"- and the one 
for 'Chemo therapy' is jF^^'-"^-measurable; the latter can be seen, using (5-separation, by 
checking that 'Chemo therapy' is locally independent from A ='Surgery' given {'Chemo 
therapy', 'Toxic reaction' } (the relevant moral graph happens to be the same as Figure 

Elb)). 
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3.4 Extensions 



Local independence graphs can easily be extended to include time-fixed covariates, such 
as sex, age and socioeconomic background of patients etc. The filtration at the start, 
JFq, then has to be enlarged to include the information on these variables. They can be 
represented by additional nodes in the graph with the restrictions that the subgraph 
on the non-dynamic nodes must be a DAG (or chain graph, cf. Gottard, 2002) and no 
directed edges are allowed to point from processes to time-fixed covariates. A process 
being locally independent of a time-fixed covariate means that the intensity does not 
depend on this particular covariate given all the other covariates and information on the 
past of all processes. (5-separation can still be applied to find further local independencies. 

The nodes in a local independence graph do not necessarily have to stand for only 
one mark (or the associated counting process), marks can be combined into one node as 
has been done in Examples 2 and 3 with 'Health status' and 'Tumour size'. This might 
be of interest when there are logical dependencies. For example if a particular illness 
is considered then the events 'falling ill' and 'recovering' from this illness are trivially 
locally dependent. If one node is used to represent a collection of marks corresponding to 
a multivariate subprocess of the whole multivariate counting process then an arrow into 
this node will mean that the intensity of at least one (but not necessarily all) of these 
marks depends on the origin of the arrow. However, some interesting information could 
be lost. For instance, if events such as 'giving birth to 1st child', 'giving birth to 2nd 
child' etc. are considered, it might be relevant whether a woman is married or not when 
considering the event of 'giving birth to 1st child' but it might not be relevant anymore 
when considering 'giving birth to 2nd child.' 

In many data situations the mark space is not finite, e.g. when measuring the mag- 
nitude of electrical impulse, the amount of income in a new job or the dosage of a drug. 
One could then discretise the mark space e.g. in 'finding a well paid job' and 'finding a 
badly paid job'. However, it must be suspected that too many of these types of events will 
generate too many logical dependencies that are of no interest and will make the graphs 
crammed. 

As we have seen in some of the examples it is sometimes sensible to consider stopped 
processes in order to avoid having to represent logical and uninteresting dependencies. 
More generally one might want to relax in Definition 12.31 the requirement 'for all t G T' 
and instead consider suitably defined intervals based on stopping times. For example, 
it might be the case that the independence structure is very different between the time 
of finishing education and starting the first job than before or after that. This deserves 
further investigation. 
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4 Discussion and conclusions 



The main point of graphical models is that they allow certain algebraic manipulations 
to be replaced by graphical ones. In the case of local independence graphs, we can read 
properties of intensity processes with respect to different, in particular reduced filtrations, 
off the graph without the need to derive explicit formulae for these intensities, and simi- 
larly we can read off relations among subprocesses such as properties (fTTl) and f|T2|) . This 
facilitates reasoning about complex dependencies, especially in the face of unobservable 
information, and simplifies calculations by reducing dimensionality. 

Clearly, it is tempting to interpret local independence graphs causally. However, we 
regard causal inference as a topic of its own and it is not the aim of this paper to go into 
much detail in this respect, but for the following few comments. Local independence graphs 
represent (in) dependencies in E{Nk{dt)\J^t-), where conditioning is on having observed J^t~ 
and, as we saw in Example 2, it makes a difference to what dependencies there are whether 
we condition on JF^- or different subsets (or even extensions) thereof. Causal inference is 
about predicting Nk{dt) after intervening in JF^-, e.g. by modulating the times of the home 
visits to be once a week in the home visits example. It is well known that conditioning on 
observation is not the same as conditioning on intervention ( "seeing" and "doing" in Pearl 

(2000) ). Hence, without further assumptions, the arrows in local independence graphs do 
not necessarily represent causal dependencies — the intensity of an event being dependent 
on whether another event has been observed before does not imply causation in the same 
way as correlation does not imply causation. Such further assumptions could be that 
all 'relevant' events (or processes) have been taken into account, like originally proposed 
by Granger in order to justify the use of the term 'causality' for what is now known as 
Granger-causality. E.g. if, in Figure [2](b), we are satisfied that by including 'Health' all 
relevant processes have been taken into account, then we could say that homevisits are 
indirectly causal for 'Death'. Obviously, in this particular example, there are many other 
relevant processes, like the occurrence of illnesses or death of the partner, that might 
be relevant. However, the literature on (non-dynamic) graphical models and causality 
has shown that causal inference is possible under weaker assumptions. Analogous results 
based on local independence graphs would require more prerequisites than we have given 
in this paper. Hence this is a topic for further research. For non-graphical approaches 
to causal reasoning in a continuous time event history setting confer Eerola (1994), Lok 

(2001) , Arjas and Parner (2004) and Fosen et al. (2004). 

Another issue is the question of statistical inference for local independence graphs. 
This can be subdivided into (i) inference when a graphical structure is given, e.g. from 
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background knowledge, but we still want to quantify the strength of the dependencies, 
and (ii) finding the graph from data if nothing about the local independence structure is 
known beforehand, which can be regarded as a particular kind of model selection or search 
task. The former has partly been addressed in section 13.31 where more specific results will 
depend on the actual modelling assumptions about the intensity processes which in turn 
will depend on the particular application. Estimation and testing within the class of 
Markov processes is tackled in Didelez (2007). More generally, local independence graphs 
can be combined with non-, semi- or parametric methods but more research is required 
to investigate how the graphical representation of the local independence structure can 
simplify inference in particular settings. As to model search, Nodelman et al. (2003) pro- 
vide a first attempt, restricted to Markov processes, at exploiting the graphical structure 
to find the graph itself when it is not postulated based on background knowledge. Clearly, 
generalisations would be desirable. 

A Appendix 

The appendix is targeted at proving Theorem 13. 5[ but to do so we first give some more 
results on the properties of local independence and (5-separation which will be used in 
that proof. These are explored along the lines of the graphoid axioms (Dawid, 1979; Pearl 
and Paz, 1987; Pearl, 1988; Dawid, 1998) which have been generalised to the asymmetric 
case by Didelez (2006). 

A.l Properties of local independence 

Proposition A.l Properties of local independence 
The following properties hold for local independence: 

(i) left redundancy: for all A, B G V: A -/^B\A, 

(ii) left decomposition: for all A,B,C CV and D C A: if A -/^B\C then D -hB\C, 

(iii) left weak union: for all A, 5, C C 1/ and C A: if A ^B\C then A^B\{CUD) 
and 

right weak union: for all A, 5, C C V" and D C B: if A 7^5|C then A -/^B\{CUD), 

(iv) left contraction: 

for all A,B,C,Dc V: ifA-f^B\C and D ^B\{AUC) then (AUD) -hB\C, 

(v) right intersection: 

for all A, 5, C C V: if A i^B\C and A i^C\B then A ^{B VJ C)\{B n C). 
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Proof: 

(i) Left redundancy holds since obviously the jF^'^'-'-^-intensities of are jF^'^^^-measurable, 
i.e. if the past of is known, then the past of is of course irrelevant. 

(ii) Left decomposition holds since the jF^'^'-'^^'^-intensities Xk{t), k ^ B, are J^^^^- 
measurable by assumption so that the same must hold for the "^-intensities Xk{t), 
k eBJoT D C A. 

(iii) Left and right weak union also trivially hold since adding information on the past of 
components that are already uninformative (left) or included (right) does not change the 
intensity. 

(iv) Left contraction holds since we have that the jFj'^^'^'-"^'-'-^-intensities A^, k G B, are by 
assumption jF/^'-'-^^'-^-measurable and these are again by assumption jF^^^'^-measurable. 

(v) The property of right intersection can be checked by noting that in the definition of 
local independence the filtration w.r.t. which the intensity process should be measurable 
is always generated at least by the process itself. □ 



Note that left redundancy, left decomposition and left contraction imply that 

A^B\C ^ A\C^B\C. (13) 
It is also always true that A -f^B \ C A -/^B\C\C, but we do not have equivalence here. 



The following property will be important for the equivalence of pairwise, local and 
global dynamic Markov properties (just like in the well-known case of undirected condi- 
tional independence graphs, cf. Lauritzen, 1996). 

Proposition A. 2 Left intersection for local independence 

Under the assumption of local independence satisfies the following property called left 
intersection: for all A, B,C G V 

if A I C and C I A then {A UC)-f^B\{An C). 

Proof: 

Left intersection assumes that the jF/^^^^'^-intensities \k{t), k G B, are J-'j^^^- as well as 
jF^'^'-'-^-measurable. With ^ we get that they are jFj^^'-^^'^-'-measurable which yields the 
desired result. 



The following can be regarded as an alternative version of the above property of left 
intersection. With f[T5]l . left decomposition and left intersection we have that for disjoint 
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A,B,C,DcV 

A -j-^B \{CVJD) and C -j^B \ {A\^ D) ^ {A\^ C) -f^B \ D (14) 
This follows from Corollary 4.3 of Didelez (2006). 

The last property that I want to consider, the 'right' counterpart of left decomposition 
given above, makes a statement about the irrelevance of a process after discarding part 
of the possibly relevant information N^^^- If the irrelevance of N^i is due to knowing the 
past of N^\/5 then it will not necessarily be irrelevant anymore if the latter is discarded. 

Proposition A. 3 Conditions for right decomposition of local independence 
Consider a marked point process and assume that the cumulative counting process ^ is 
non-explosive and that intensites exist. Let A, B,CgV,DgB, with {Br\A)\{CUD) = 
0. The following property, called right decomposition, 

A-f^B\C A -f^D I C 

holds under the conditions that 

B i^A\{CUD) I (CUD) 

and 

A ^{k} \CUB or B ^{k} \{CUDUA) 

for all k e C\D. 
Proof: 

In this proof we proceed somewhat informally for the sake of simplicity. The formal proof 
is based on the results of Arjas et al. (1992) and given in Didelez (2000, pp. 72). 

Redefine A* = A\{C U D), B* = B\D and C* = C\D. Then A* n B* = (D and with 
(fT2|) and (fTTl) . it can be shown that the assumptions of the present proposition imply that 

as well as 

J^/^*Xjf * I Jf^^. (16) 
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We want to show that the J^/^^'^^'^-intensity Ad(^) of N£,(t) is ^-^-measurable. With 
the above and interpretation ([T]) we have 

~XD{t)dt = E{ND{dt) I = E{ND{dt) \ J^fy^^^^) 

= E{E{Nn{dt) I jrA*uB'uC'uD^ | ^a*uc*uD) 
= E{E{ND{dt) I J^fS^c'uD^ I jrA*uc'uD^ ^gj^g ^ 
= E{ND{dt) I J^,*!*^^) using ([ISD 

= E{Nn{dt)\J^^y''), 

as desired. 

A. 2 Properties of (5— separation 

For a general investigation of the properties of (5-separation we need to complete Defini- 
tion [XSl by the case that A,B and C are not disjoint: we then define that C (5-separates 
A from B if C\B (5-separates A\{B U C) from B. We further define that the empty set 
is always 5-separated from B. Additionally, we define that the empty set ^-separates A 
from B a A and B are unconnected in (G'An(AuB))™' 

It can be shown (Didelez, 2006) that 5-separation satisfies the same properties as lo- 
cal independence given above in Proposition lA. ll if we replace A -/^B\C hj "C (5-separates 
A from i?" which we will write as AlRsB\C. In particular it satisfies left redundancy, left 
decomposition, let and right weak union, left and right contraction as well as left and 
right intersection without requiring further asumptions. The property of right decompo- 
sition holds for ^-separation under conditions analogous to those in Proposition IA.3I In 
particular we have the following special case of right decomposition 

AiRsB \C, D c B ^ AiRsD \ {CU B)\D (17) 

which is Lemma 4.11 in Didelez (2006). 

In addition, we want to show how (5-separation can be read off a local independence 
graph in a different but equivalent way to Definition 13.31 We mention this firstly, because 
it will be more familiar to readers who use ^-separation for DAGs (Pearl, 1988; Verma 
and Pearl, 1990), and secondly, because some parts of the proof of Theorem 13.51 are easier 
to show using this alternative way of checking (5-separation. 

First, the different notions of paths and trails need to be made more stringent. Con- 
sider a directed or undirected graph G = {V, E). An ordered {n + l)-tuple (jo, • • • , jn) of 
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distinct vertices is called an undirected path from jo to j„ if {ji-i,ji} G E and a directed 
path if jj) G for all i = 1, . . . ,n. A (directed) path of length n with jo = jn 
is called a (directed) cycle. A subgraph vr = of G with y = {jo, •••,jn} and 

E' = {ci, . . . , e„} C -E is called a traz/ between jo and j„ if = (ji-iji) or = (ji, jj-i) 
or Cj = {ji, jj-i} for all z = 1, . . . , n. Further, for a directed graph we say that a trail 
between j and is blocked by C if it contains a vertex 7 such that either (i) directed 
edges of the trail do not meet head-to-head at 7 and 7 G C, or (ii) directed edges of the 
trail meet head-to-head at 7 and 7 as well as all its descendants are no elements of C. 
Otherwise the trail is called active. 

Proposition A. 4 Trail condition for 6-separation 

Let G = {V, E) be a directed graph and A, B, G pairwise disjoint subsets of V. Define 
that any allowed trail from A to B contains no edge of the form {b,k),b & B,k E V\B. 
For disjoint subsets A, B, G of V, we have that G 5-separates A from B if and only if all 
allowed trails from A to B are blocked by G. 

The proof is given in Didelez (2000, pp.22; cf. also Didelez, 2006). 
A.3 Proof of Theorem [331 

It is easily checked that ([7]) ^ ([5]) ^ ([1]): First, pa(/c) always 5-separates V\{pa,{k) U {k}) 
from {k} in G, hence ([5]) is just a special case of ([7]). Also, it is easy to see that the 
equivalence of the pairwise and local dynamic Markov properties immediately follows from 
left intersection assuming left weak union and left decomposition. Thus, the following 
proof considers situations where A, B, G do not form a partition of V or pa(i?) (/ G. The 
structure of the proof corresponds to the one given by Lauritzen (1996, p. 34) for the 
equivalence of the Markov properties in undirected conditional independence graphs. Due 
to the asymmetry of local independence, however, this version is more involved. 
Assume that (jlj) holds and that G 5-separates A from B in the local independence graph. 
We have to show that A ^ B \ G, i.e. the J^^^^-^^'^-intenities Afc(t), k e B, are J^f 
measurable. The proof is via backward induction on the number \G\ of vertices in the 
separating set. If \G\ = \V\ — 2 then both, A and B consist of only one element and ([7]) 
trivially holds. If |C| < |y| — 2 then either A or B consist of more than one element. 
Let us first consider the case that A,B,G is a partition of V and none of them is empty. 
If > 1 let a G A. Then, by left weak union and left decomposition of 5-separation we 
have that G U (A\{a}) 5-separates {a} from B, i.e. 

{ajlRsB \ GU{A\{a}) 
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and C U {a} 5-separates from B in G, i.e. 

A\{a}msB I (CU {a}). 

Therefore, we have by the induction hypothesis that 

{a} -hB\CU {A\{a}) and A\{a} -f^B\{CU {a}). 

From this it follows with the modified version of left intersection as given in f[T^ (which 
can be applied because of the assumption that ([6]) holds) that A ^B \C a.^ desired. 
If \B\ > 1 we can show by a similar reasoning, applying (fT7j) to {(3} G -B, that A -j^B \ C . 
Let us now consider the case that A^B^C C V are disjoint but no partition of V . First, 
we assume that they are a partition of An(74 U i? U C), i.e. that A U i? U C is an ancestral 
set. Let 7 G V^\(^ U -B U C), i.e. 7 is not an ancestor oi AVJ B VJ C . Thus, every allowed 
trail (cf. Proposition IA.4P from 7 to i? is blocked hy AVJ C since any such trail includes 
an edge (/c, h) for some h E B where no edges meet head-to-head in k and /c G A U C 
Therefore, we get 

{^]msB I (AUC). 

Application of left contraction, weak union, and decomposition for (5-separation yields 

AlRsB I (CU{7}). 

It follows with the induction hypothesis that 

A-f^B\{CU {7}) as well as {7} -f^B\{AU C). 

With left intersection as given by fll4p and left decomposition for local independence we 
get the desired result. 

Finally, let A, B, C be disjoint subsets of V and AVJ B U C not necessarily an ancestral 
set. Choose 7 G an(A UBUC) and define = GAn(AuBuc)- ^^^^^ B\Cm (G^)"* 
we know from the properties of ordinary graph separation that 

(i) either {-f}ALg B\{AuC)m {G^)"" 
(li) or AALy {7} {(BUG) in (G^)™. 

(i) In the first case {7} iRsB \ {AU G) and it follows from left contraction that 

{A U {7}) iRsB I G. 

Application of left weak union and left decomposition yields AlRsB \ {GU {7}). With the 
induction hypothesis we therefore get 

A^B\{GU {7}) and {7} ^B \ {AU G). 
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Left intersection according to (fT^ and left decomposition for local independence yields 
A-f^B\C. 

(ii) The second case is the most complicated and the proof makes now use of right de- 
composition for local independence under the conditions given in Proposition IA.3I First, 
we have from (ii) that y4_LLc,{7} | B U C m. (G'An(Aufiuc))'" since the additional edges 
starting in B can only yield additional paths between A and 7 that must be intersected 
by B. Since deleting further edges out of 7 does not create new paths, it holds 

AiRsii} I (BUC). 

With AlRsB I C, application of right contraction for 5-separation yields A\Rs{B U {7}) | 
C. Now, we can apply property (fT7|) to get AlKsB \ {C U {7}) from where it follows with 
the induction hypothesis that 

A 7^5 I (C U {7}) and A 7^{7} | {B U C). 

With right intersection for local independence we get A -J^ {B VJ {7}) | C. In addition, 
{7} A \ (i? U C) by the same arguments as given above for A -j^ {7} | (5 U C). In 
order to apply Proposition IA.3I we still have to show that for all A; G C either A\R^{k\ \ 
[C \^ B \^ {7}) or {7} iR^jA;} | (C U -B U A) which by the induction hypothesis implies 
the corresponding local independencies. To see this, assume that there exists a vertex 
/c G C for which neither holds. With the trail condition we then have that in G'An(AuBuc) 
there exists an allowed trail from A and 7, respectively, to k such that every vertex where 
edges do not meet head-to-head are not in (C U i? U {7})\{/i;} and (C U -B U 
respectively, and every vertex where edges meet head-to-head or some of their descendants 
are in {CV^B\^ {7})\{/c} respective (C U 5 U A^)\{}^\. This would yield a path between A 
and 7 which is not blocked by C U 5 (note that A; is a head-to-head node on this trail) in 
^An(yiuBu(7) • This in turn contradicts the separation of A and 7 by i? U C in G^p^^i^j^^j^yjc) 
because the edges starting in B cannot contribute to this trail. Consequently we can apply 
right decomposition and get the desired result. 

A.4 Proof of Theorem 13.61 

The factorisation ffTOl) implies that the marginal likelihood for the marked point process 
discarding events not in An(y4 U i? U C), i.e. {(T^, | s = 1, 2, . . . ; G i^An(AuBuc)}5 is 
given by 

Lit I i/AnCAuSuC)^ ^ -Q ^^^^ I ^cl(fc)) 

feeAn(AuBuC) 
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as none of the intensities of A^^, k ^An{A U B U C) depend on V\An{A U B U C). Hence, 
the hkehhood may be written as product over factors that only depend on cl{k), k G 
An{A U BUC). Let C = {d{k) \ k e An{A U 5 U C)} be the set containing all such sets 
and let gc{t\-), c G C, be these factors. Then, we have 

Further, by rearranging the factors the sets in C can be taken to be the 'cliques', i.e. 
the maximal fully connected sets of nodes, of the graph (GAn(AuBuc))'"- The above thus 
corresponds to the factorisation property of undirected graphs which in turn implies the 
global Markov property for undirected graphs (Lauritzen, 1996, p. 35). This means that 
when we have a separation, like C separating A and B, in this graph (G'An(AuBuc))™ ^^e 
corresponding conditional independence fill I) holds, which completes the proof. 



Acknowledgements: 

I wish to thank Iris Pigeot for her persistent encouragement and advice, and Niels Keiding 
for valuable comments and suggestions. Financial support by the German Research Foun- 
dation (SFB 386) is also gratefully acknowledged. The referees' comments have greatly 
helped to improve the paper. 

References 

Aalen, O. O. (1987) Dynamic Modelling and Causality. Scandinavian Actuarial Jour- 
nal, 177 - 190. 

Aalen, O. O., 0. Borgan, N. Keiding, and J. Thormann (1980) Interaction between Life 
History Events: Nonparametric Analysis of Prospective and Retrospective Data in 
the Presence of Censoring. Scandinavian Journal of Statistics 7, 161 - 171. 

Allard, D., Brix A. and J. Chadoeuf (2001) Testing local independence between two 
point processes. Biometrics 57 508 - 517. 

Andersen, P. K., 0. Borgan, R. D. Gill, and N. Keiding (1993) Statistical Models Based 
on Counting Processes. Springer, New York. 

Arjas, E., P. Haara, and I. Norros (1992) Filtering the Histories of a Partially Observed 
Marked Point Process. Stochastic Processes and Applications 40, 225 - 250. 

Arjas, E., and J. Parner (2004) Causal Reasoning from Longitudinal Data. Scandinavian 
Journal of Statistics 31, 171 - 201. 



27 



Bremaud, P. (1981) Point Processes and Queues: Martingale Dynamics. Springer, New 
York. 

Comte, F., E. Renault (1996) Long memory continuous time models. Journal of Econo- 
metrics 73, 101 - 149. 

Cowell, R. G., A. P. Dawid, S. L. Lauritzen, and D. J. Spiegelhalter (1999) Probabilistic 
Networks and Expert Systems. Springer, New York. 

Cox, D. R. and V.S. Isham (1980) Point Processes. Chapman and Hall, London. 

Cox, D. R. and N. Wermuth (1996) Multivariate Depencencies — Models, Analysis and 
Interpretation. Chapman and Hall, London. 

Dahlhaus, R. (2000) Graphical Interaction Models for Multivariate Time Series. 

Metnka 51, 157 - 172. 

Dawid, A. P. (1979) Conditional Independence in Statistical Theory. Journal of the 
Royal Statistical Society, Series B 4i, 1 ^ 31. 

Dawid, A. P. (1998) Conditional Independence. In S. Kotz, C. B. Read, and D. L. 
Banks (Eds.), Encyclopedia of Statistical Sciences, Update Volume 2, pp. 146 - 155. 
Wiley-Interscience, New York. 

Dean, T. and K. Kanazawa (1989)) A model for reasoning about persistence and cau- 
sation. Comp. Intelligence 5, 142 - 150. 

Didelez, V. (2000) Graphical models for event history data based on local independence. 
Logos, Berlin. 

Didelez, V. (2006) Asymmetric separation for local independence graphs. In R. Dechter 
and T. Richardson (Eds.), Proceedings of the 22nd Annual Conference on Uncer- 
tainty and Artificial Intelligence, pp. 130 - 137. UAI Press, Oregon. 

Didelez, V. (2007) Graphical models for composable finite Markov processes. Scandi- 
navian Journal of Statistics 34, 169 - 185. 

Edwards, D. (2000) Introduction to Graphical Modelling (2nd ed.) Springer, New York. 

Eerola, M. (1994) Probabilistic Gausality in Longitudinal Studies, Volume 92 of Lecture 
Notes in Statistics. Springer, Berlin. 

Eichler, M. (1999) Graphical Models in Time Series Analysis. Ph. D. thesis. University 
of Heidelberg. 

Eichler, M. (2000) Causality Graphs for Stationary Time Series. Technical report. Uni- 
versity of Heidelberg. 



28 



Fleming, T. R. and D. P. Harrington (1991) Counting Processes and Survival Analysis. 
Wiley, New York. 

Florens, J. P. and D. Fougere (1996) Noncausality in Continuous Time. Economet- 
rica 61 1195 - 1212. 

Florens, J. P., Mouchart, M. and J. M. Rolin (1990) Elements of Bayesian Statisitcs. 
Dekker, New York. 

Fosen, J., Borgan, O., Weedon-Fekjoer, H. and O. O. Aalen (2004) Path analysis for 
survival data with recurrent events. Technical Report, University of Oslo, Norway. 

Granger, C. W. J. (1969) Investigating Causal Relations by Econometric Models and 
Cross-spectral Methods. Econometrica 37, 424 - 438. 

Gottard, A. (2007) On the inclusion of bivariate marked point processes in graphical 
models. Metrica 66, 269 - 287. 

Jordan, M. (1999) Learning in Graphical Models. MIT Press, Cambridge. 

Keiding, N. (1999) Event history analysis and inference from observational epidemiol- 
ogy. Statistics in Medicine 18, 2353 - 2363. 

Koster, J.T.A. (1996) Markov properties of non-recursive causal graphs. Annals of 
Statistics 24, 2148 - 2177. 

Lauritzen, S. L. (1996) Graphical Models. Clarendon Press, Oxford. 

Lok, J. J (2001) Statistical modelling of causal effects in time. Ph. D. thesis. Free 
University of Amsterdam. 

Nodelman, U., Shelton, C. R. and D. Koller (2002) Continuous Time Bayesian Net- 
works. Proceedings of the 18th Annual Conference on Uncertainty in Artificial In- 
telligence, 378 - 387. Morgan Kaufman, San Francisco. 

Nodelman, U., Shelton, C. R. and D. Koller (2003) Learning Continuous Time Bayesian 
Networks. Proceedings of the 19th Annual Conference on Uncertainty in Artificial 
Intelligence, 451 - 458. Morgan Kaufman, San Francisco. 

Pearl, J. (1988) Probabilistic Reasoning in Intelligent Systems. Morgan Kaufman, San 
Mateo. 

Pearl, J. (2000) Causality. Models, Reasoning, and Inference. Cambridge University 
Press, Cambridge. 

Pearl, J. and A. Paz (1987) Graphoids: A Graph-Based Logic for Reasoning about 
Relevancy Relations. In B. Du Boulay, D. Hogg, and L. Steel (Eds.), Advances in 
Artificial Intelligence - //, pp. 357 - 363. 



29 



Robins J. (1986) A new approach to causal inference in mortality studies with sustained 
exposure period — applications to control of the healthy workers survivor effect. 
Mathematical Modelling 7, 1393 - 1512. 

Robins J. (1997) Causal inference from complex longitudinal data. In M. Berkane, 
Latent Variable Modeling and Applications to Causality. Lecture Notes in Statistics 
120. Springer, New York. 

Schweder, T. (1970) Composable Markov Processes. Journal of Applied Probability 7, 
400 - 410. 

Spirtes, P. (1995) Directed cyclic graphical representations of feedback models. In Pro- 
ceedings of the 11th Annual Conference on Uncertainty and Artificial Intelligence, 
pp. 491 - 498. Morgan Kaufman, San Francisco, California. 

Spirtes, P., C. Glymour, and R. Scheines (2000) Causation, Prediction, and Search (2nd 
ed.) MIT Press, Cambridge, Massachusetts. 

Vass, M., Avlund, K., Hendriksen, C, Andersen C. K. and N. Keiding (2002) Preventive 
home visits to older people in Denmark: methodology of a randomized controlled 
study. Aging Clin. Exp. Res 14, 509 - 515. 

Vass, M., Avlund, K., Kvist, K., Hendriksen, C, Andersen C. K. and N. Keiding (2004) 
Structured home visits to older people. Are they only of benefit for women? A 
randomised controlled trial. Scand. J. Prim. Health Care 22, 106 - 111. 

Verma, T. and J. Pearl (1990) Causal Networks: Semantics and Expressiveness. In 
R. D. Shachter, T. S. Levitt, L. N. Kanal, and J. F. Lemmer (Eds.), Uncertainty 
and Artificial Intelligence, pp. 69 - 76. North Holland, Amsterdam. 

Whittaker, J. (1990) Graphical Models in Applied Multivariate Statistics. Wiley, Chich- 
ester. 



30 



